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Abstract 


The  Weibull  distribution  model  for  brittle  fracture 
is  applied  to  a  hole  size  study  in  uniaxial  tension  and 
to  a  study  in  biaxial  tension.  This  application  of  the 
Weibull  model  uses  a  numerical  integration  of  stress 
functions  across  a  high-failure-risk  volume.  The  results 
for  uniaxial  tension  indicate  a  hole  size  effect  which 
agrees  in  form  with  the  "point-stress"  or  characteristic 
dimension  theory .  The  uniaxial  strength  predictions 
based  on  Weibull  theory  are  uniformly  conservative.  The 
predictions  also  show  increasing  notch  sensitivity  as 
the  material  becomes  more  "perfect",  that  is,  has  fewer 
and  smaller  inherent  flaws.  The  biaxial  Weibull  study 
accurately  predicts  the  failure  mode  and  the  strengthen¬ 
ing  effect  of  biaxial  tension.  The  biaxial  strength  pre¬ 
dictions  are  also  generally  conservative. 
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Author's  Note  on  Weibull  Statistics 

It  is  important  for  the  reader  to  realize  the  difference 
between  the  use  of  Weibull  statistics  for  curve  fitting  and 
the  calculation  of  "Weibull-type  integrals." 

The  Weibull  distribution  may  be  used  purely  as  a  curve- 
fit  scheme.  In  this  use,  it  has  general  utility,  much  as  a 
normal  distribution  or  any  other  distribution  for  describing 
data  scatter.  This  paper,  however,  uses  a  much  more  specific 
application,  that  of  "Weibull-type  integrals."  In  this  appli¬ 
cation,  the  parameters  of  the  Weibull  distribution  have 
physical  significance;  they  are  related  to  stress  and  geometry 
variables.  These  integrals  attempt  to  calculate  parameters 
for  a  Weibull  distribution,  rather  than  simply  infer  them 
from  data  as  in  a  curve-fit  scheme. 
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Introduction 


The  problem  of  strength  degradation  in  the  presence 
of  a  notch  is  important  to  the  materials  designer.  These 
notches  may  be  introduced  intentionally,  as  in  the  case 
of  bolt  holes  or  cutouts,  or  unintentionally,  as  in  the 
case  of  a  materials  fabrication  error  or  a  fatigue  crack. 

The  designer  needs  formulae  which  help  him  estimate 
strength  "knock-down  factors"  for  his  design;  he  also 
needs  a  way  of  comparing  the  notch  sensitivity  of  dif¬ 
ferent  materials.  Predictions  of  notch  sensitivity  at 
a  variety  of  hole  sizes  based  on  the  finite  state  of 
stress  in  the  material  have  been  successful, 

iPipes,  Wetherhold,  Gillespie] , [Pipes ,  Gillespie,  Wetherhold] . 
This  model  is  known  also  as  the  "point-stress"  model. 

The  meaning  of  the  notch  sensitivity  parameters  found 
in  the  above  references  has  not,  however,  been  related 
to  any  specific  physical  variables. 

Another  important  factor  for  a  materials  designer 
is  to  understand  how  the  inherent  variability  of  his 
material  will  affect  design.  Factors  of  reliability  must 
often  be  calculated  based  on  a  number  of  strength  data 
to  insure  a  certain  low  probability  of  failure.  Weibull 
statistics  are  often  employed  as  a  good  way  of  calculating 
reliability. 
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This  paper  joins  the  strength  degradation  predictions 
of  the  point-stress  model  to  the  use  of  Weibull  statistics . 
In  this  way,  the  parameters  contained  in  the  point-stress 
strength  equations  are  related  to  the  unnotched  tensile 
strength  data  scatter.  Thus,  the  variability  in  unnotched 
strength  is  related  to  the  strength  degradation  in  the 
presence  of  a  notch  or  flaw.  The  specific  calculations 
in  this  paper  are  for  both  through-holes  and  slits  (notches) 
in  a  uniaxial  tension  plate,  and  for  through-holes  in  a 
biaxial  tension  plate. 


iii 


Part  1  UNIAXIAL  TENSION 


1.1  Weibull  Theory 

The  Weibull  brittle  failure  model  may  be  used  to 
predict  conservative  failure  loads  for  notched  composites. 
The  prediction  is  based  on  an  integration  which,  ideally, 
would  be  over  the  entire  stressed  volume.  For  a  plane 
stress  state  of  a  notched  plate  with  one  predominant 
stress  direction  "x,"  the  stress  state  may  be  given  by 

<^jj(x,y)  =  t(x,y)  1.1 

where  a”  is  the  far  field  stress,  and  t(x,y)  is  a 
geometric  locus,  [see  Appendix  B] .  See  Fig.  la.  In  general, 
t(x,y)  would  come  from  a  complex  variable  solution  for 
stress.  Fig.  la’  shows  an  elliptic  center  through-notch, 
although  the  actual  notch  could  be  any  shape . 

Using  the  results  of  Appendix  A,  we  note  that  the 
average  far-field  failure  stress  can  be  found  as  follows. 

The  stress  integral  used  in  reliability  calculations  is 

B(a“)  =  /  dV  1.2 

X  ^6 

where  we  assume  that  a  =  a  is  the  predominant  stress. 

Since  equation  1.1  introduces  the  geometric  variation  in 
stress,  we  may  write 
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1.3 


Note  that  a”  may  be  removed  from  the  integral.  In  order 
to  place  equation  1.3  in  a  standard  form,  we  define  a 
new  notch  scale  parameter  (dependent  on  notch  geometry) 
by 

\  -1/a 

/  [t(x,y)°^]dV  )  6  1.4 

V  / 


This  way,  the  reliability  function  R  for  probability  of 
survival  at  any  given  far-field  stress  c“  is 

1.5 

R(a°°)  =  ^  1.6 


R(a  )  =  exp 


N 


a 


The  average  value  of  far-field  stress,  the  most 
meaningful  statistic  from  small  quantity  notched  strength 
tests,  can  now  be  calculated.  From  Appendix  A. I,  the 
average  notched  strength  is  related  to  the  scale 
parameter  by 

Fjj  =  Bjj  r  (1+1/a)  1.7 


We  now  hypothesize  that  the  volume  of  integration  can 
be  reduced  to  a  thin  strip  perpendicular  to  the  applied 
loading.  See  Fig.  lb.  The  strip  has  width  "6”  by  plate 
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thickness  "h. 


"  Within  the  strip  “-j  £  Y  , 

t(x,y)  =  t(x,0)  (stress  within  the  volume  =  stress  at  the 
line  Y  =  0) ;  6  £  .05  a  is  sufficiently  small  for  error  to 
be  _<  1%.  The  strip  length  could  be  from  Y  =  a  to  the  plate 
edge  (finite  plate) ,  or  to  10  times  maximum  a  (infinite 
plate) .  The  result  of  this  narrowing  of  the  integration 
volume  of  interest  is  to  reduce  equation  1.4  to 

=  (  h6  ^  [t(0,y)]“  dY)  3  1.8 

We  would  expect  such  a  calculation  to  give  conservative 
results,  as  the  volume  under  consideration  is  the  highest 
risk  vlume.  The  average  predicted  failure  load  can  be 
found  by  combining  equation  1 . 8  with  1.7, 


In  order  to  calculate  the  notched/unnotched  strength 
ratio,  we  make  the  following  assumption; 

a,  the  shape  parameter, is  identical 
for  the  notched  and  unnotched  specimens. 

The  average  strength  for  an  unnotched  specimen,  using 
the  same  volume  of  integration  as  for  the  notched  specimen 
(with  a  =  0)  is 

=  [h6L]"^/“  3r (1+1/a)  1.10 
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This  leads  to  an  equation  predicting  the  reduction 
in  strength  with  a  notch;  dividing  equation  1.9  by  1.10, 


/  ^ 

a  \ 

[t(0,Y)]°‘dy 

^N  ^  ^ 

/ 

L-l/a 

1.11 

where  we  omit  the  super-bar,  but  define  the  lefthand  side 
as  average  values.  Note  that  as  a  0,  ->  1  as  it 

should.  We  are  now  in  a  position  to  vary  the  parameters 
in  equation  1.11,  evaluating  the  integral  1.11  numerically. 

1.2  Point  Stress  Theory 

This  failure  theory  is  introduced  by  manner  of 
comparison  with  Weibull  theory.  The  size  effect  predicted 
by  point  stress  theory  is,  in  fact,  identical  in  form  to 
that  of  Weibull  theory. 

The  loading  geometry  is  identical  to  Fig.  1.  The 
point  stress  failure  prediction  is  that  when  the  stress 
at  some  distance  "<1^"  into  the  material  reaches  the  unnotched 
ultimate  stress,  the  sample  fails,  i.e.  failure  occurs  when 


where  is  the  unnotched  ultimate  [Whitney,  Nuismer] .  To 
provide  an  absolute  hole  size  effect  in  agreement  with 
experiment  to  vary  with  notch  size  per  [Pipes,  et.al.] 


5 


1.13 


=  (a/a^)"‘/K 

where  a  =  notch  half-length 

^  =  one  times  units  used  (arbitrary) 

K  =  notch  sensitivity  parameter 
m  =  exponential  parameter 
Since  the  direction  of  fracture  is  presumed  perpendicular 
to  the  applied  load,  only  the  stress  formula  for  a(x=0,y) 
is  required.  The  reference  notch  half  size  ^  is  included 
only  to  avoid  awkwardness  in  the  units  of  K. 
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1.3  Numerical  Integral  Results 


The  stress  profile  in  the  selected  volume  of  integration 
is  complex  enough  that  equation  1.11  demands  ntamerical 
integration.  The  integration  routine  selected  is  DCADRE, 
a  part  of  the  International  Mathematical  and  Scientific 
Library  (IMSL) . 

1.3.1  Circular  Hole 

For  a  circular  through  hole  geometry,  we  have  the 
infinite  plate  stress  profile  [Konish,  Whitney] 

t(0,Y)  =  1+|p^-8-|p^-(k”-3)  (|p^-^^)  1.14 

n 

where  p  =  ^ 

R  =  hole  radius 

k”  =  classical  stress  concentration  factor, 

calculatable  from  the  elastic  constants 
R  1 

For  a  finite  plate  where  ^  <  j,  we  may  use  an  isotropic 
finite  width  correction  (FWC)  and  multiply  the  results  of 
equation  1.14  by  the  FWC  factor  [Peterson] 

FWC  =1-0.05X+1.5X^  1.15 

X  =  2R/W 

W  =  sample  width 
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A  hole  size  series  was  run  in  which  the  specimen 
width  was  held  constant,  and  the  hole  radius  was  varied 
at  a  given  Weibull  shape  parameter  a.  The  material  being 
modeled  was  isotropic;  thus,  =  3. 

The  resulting  values  were  used  to  back-calculate 

a  d  value.  Since,  for  circular  notches, 
o 


N  _ 


1  ^2  .  3  ^4 


=  1  +  i  +  i  5’  -  (k:  -  3)4  5®  -  J  4) 


where 


K  = 


2  '2 
R 


1.16 


R+d 


this  implies  an  inverse  function  "g," 

fa. 


do  =  g 


;  R 


1.17 


The  values  of  d  calculated  from  the  o„/a  values  were 
o  N  o 

then  fitted  to  the  equation  1.13,  taking  logarithms: 


In  do  =  m  In  (:^)  -  In  k  1.18 

^o 

The  results,  shown  in  Table  1,  demonstrate  that  the 
presumed  form  for  d^  as  a  function  of  radius  (1.13  or 
1.18)  is  extremely  accurate. 


Increasing  a  is  also  seen  to  bring  increased  notch 
sensitivity.  Both  k  and  m  increase  with  increasing  a,  which 
indicates  a  more  rapid  strength  degradation  in  the  presence  of  a 
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notch,  (see  Table  1,  Fig.  2,3)  Physically,  a  material 
with  high  a  has  a  very  tight  distribution  of  small  flaws . 
Thus,  when  we  add  a  notch,  it  far  outweighs  any  innate 
flaw  effects.  Conversely  for  low  a  there  are  various 
size  flaws  in  the  material,  and  much  data  scatter.  The 
addition  of  one  more  flaw  (the  notch)  is  almost  unnoticed; 
notch  sensitivity  (m  and  k)  are  low.  Note  that  the  curves 
in  Fig.  2,3  have  a  "kneeing  over";  notch  sensitivity 
reaches  a  rough  constant  for  a  greater  than  about  30. 


Table  1 

Weibull  Shape  Parameter  a  versus  notch  sensitivity 
Constants  k  and  m  for  Circular  Notch 


a 

m 

k 

r* 

3 

0.313 

3.89 

0.985 

5 

0.386 

5.49 

0.987 

10 

0.520 

9.70 

0.990 

12 

0.573 

10.75 

0.986 

15 

0.560 

14.0 

0.983 

20 

0.611 

16.8 

0.996* 

30 

0.572 

27.8 

0.977* 

*  correlation  coefficient  from  the  least  squares  fit  of 
linear  equation  „ 

In  d^  =  m  In  (^)  -  In  k; 

°  ^o 

four  points  (radii)  used  for  each  given  a. 


Same  as  *,  except  three  points  (radii)  used  for  each  given  a. 
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FIGURE  2 

m  CALCULATED  AS  A  FUNCTION  OF  a  FROM  NUMERICAL  INTEGRALS 
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CIRCULAR  HOLE  NOTCH  SENSITIVITY 


FIGURE  3 

k  CALCULATED  AS  A  FUNCTION  OF  a  FROM  NUMERICAL  INTEGRALS 


ir 


1.3.2  Slit  Notch 

For  a  through  slit  notch  (crack) ,  results  similar  to 
the  hole  results  were  obtained.  A  degradation  of  strength 
with  increasing  hole  size  is  observed,  and  the  notched 
strength  reduction  calculated  by  Weibull  integrals  is  very 
conservative . 

The  theoretical  stress  profile  [Savin]  must  be  modified 
for  use  in  the  niimerical  integral .  *  We  let 

Y  ^  I.OIC 

1.19 

1.01c  >  Y  >  C 

The  entire  profile  is  multiplied  by  a  finite  width  correction 
factor  [Paris,  Sih] 


where  W  =  sample  width 

C  =  crack  half-length 

The  integrals  of  equation  1.11  were  performed  for  a  given 
width  sample,  with  the  crack  sized  varying  at  a  given  a. 

The  numerical  integrals  were  used  to  back  calculate  a 
value  of  d^  by  inverting  the  point  stress  criterion; 

*  The  IMSL  Routine  (DCADRE)  uses  progressively  smaller  mesh 
sizes  and  tries  to  limit  the  error  between  mesh  divisions. 

The  theoretical  profile  is  infinite  at  Y=C;  while  the  integral 
exists,  this  form  is  unsuitable  for  numerical  integration. 
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since,  for  cracks. 


thus 


1.21 


1.22 


The  results  show  that  the  variation  of  d  with  crack 

o 

size. 

In  =  m  In(^)  -  In  k  1.23 

^o 

fits  the  variation  in  notch  strength  extremely  well. 

(See  Table  2) .  The  variation  in  m  and  k  with  a  may  be 
seen  in  Figures  4,  5. 


Table  2 

Weibull  Shape  Parameter  a  versus  notch  sensitivity 
Constants  k  and  m  for  Slit  Notch 


a 

m 

3 

0.189 

5 

0.448 

10 

0.723 

12 

0,758 

15 

0.801 

20 

0.836 

k 

r* 

10.9 

0.986 

23.9 

0.999 

49.3 

0.999 

56.8 

0.999 

64.7 

1.000 

75.3 

1.000 

*  Correlation  coefficient  from  the  last  squares  fit  of  the 
linear  form 

In  d^  =  m  In  (^)  -  In  k; 

Four  points  (crack  sizes)  u§ed  for  each  a. 
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FIGURE  4 

m  CALCULATED  AS  A  FUNCTION  OF  a  FROM  NUMERICAL  INTEGRALS 
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,  WEIBULL  SHAPE  PARAMETER 


001 


FIGURE . 5 

k  CALCULATED  AS  A  FUNCTION  OF  a  FROM  NUMERICAL  INTEGRALS 
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WEIBULL  SHAPE  PARAMETER 


1.4  Conclusions  and  Recommendations 

1.  The  results  of  the  Weibull  calculation  are  conservative. 
Compare  the  notch  sensitivity  constants  with  experiment 
[Whitney,  Kim  1976]  for  a  [90/0/±45]  Graphite/Epoxy: 


Notch 

Geometry 

CALCULATED 

USING  a=20 

EXPERIMENTAL 
(Curve  Fit) 

m 

k 

m 

k 

Hole 

0.61 

16.8 

0.49 

9.6 

Slit 

0.84 

75.3 

0.39 

19.4 

The  higher  values  of  k  and  m  calculated  by  Weibull 
integrals  predict  a  more  rapid  degradation  of  notched 
strength  with  larger  notch  size  than  is  actually  experienced. 
For  circular  holes,  this  conservatism  is  on  the  order  of 
a  design  safety  factor,  with  calculated  notched  strengths 
on  the  order  of  1/1.5  to  1/2  lower  than  seen  experimentally. 
For  slit  notches  the  calculation  is  extremely  conservative. 
One  explanation  for  this  is  that  the  theoretical  stress 
profile,  even  our  altered  form  1.19  is  more  strenuous 
than  that  actually  experienced.  The  model  is  an  atomically 
thin  crack;  any  machining  takes  the  form  of  a  slot. 

Other  reasons  for  conservatism  on  both  circular  and 
slit  notches  are:  The  volume  of  integration  is  the  highest 
risk  volume  (highest  stresses) .  Also,  the  Weibull  theory 
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is  a  weak  link  theory  (see  appendix  E)  where  failure  of 
any  volume  brings  about  failure  of  the  whole.  This  is  a 
correct  assumption  for  very  brittle  materials  such  as 
ceramics,  but  is  somewhat  conservative  for  composites. 

II  The  notch  sensitivity  constants  k,  m  are  volume 
dependent.  If  one  chooses  a  larger  volume,  the 
Weibull  calculation  predicts  different  notch  sensitivity 
constants  m,k  (at  a  given  a).  For  example: 


GEOMETRY 

CALC. 

a=10. 

USING 

W=2" 

CALC. 
a=10 , 

USING 
W=4  " 

m 

k 

m 

k 

Hole 

0.52 

9.7 

0.61 

6.1 

Slit 

0.72 

49.3 

0.76 

37.6 

This  suggests  that  the  sample  volume  in  a  notch-size 
series  should  be  kept  as  constant  as  possible  and  practical. 
One  suggestion  for  correcting  between  various  volumes  used 
experimentally  is  to  define  a  standard  volume,  and  scale 
experimental  results  to  the  standard  by  multiplying  by 
the  ratio  of  the  theoretical  calculations  for  the  voliimes. 
That  is,  if  we  let 

S  =  experimental,  at  volume  2 

S*  =  experimental,  "shifted"  to  standard  volume  1 
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=  c5„/a  calculated  for  volume  2  at  given  a 
2  N  o 


volume  1  at  given  a 


We  let  the  change  in  volvime  be  represented  by  a  varying 
width  "w. "  Thus, 


S*  =  S 


1.23 


Corrections  of  this  type  have  been  modest,  on  the  order 
of  a  few  percent. 


Ill  The  notch  sensitivity  of  a  material  increases  with 
increasing  a.  A  high  a  value  (from  unnotched  tensile 
tests)  implies  a  very  tight  (narrow)  distribution  of 
inherent  flaws.  The  rapid  strength  degradation  with  notch 
size  for  a  material  with  high  a  implies  that  these  inherent 
flaws  are  small.  Due  to  the  inherent  "perfection"  of  the 
material,  addition  of  a  notch  has  a  dramatic  effect  on 
strength. 

Materials  with  low  a  thus  have  lower  notch  sensitivites ; 
however,  this  is  a  mixed  blessing.  Low  a  values  imply  a 
wide  distribution  of  various  size  inherent  flaws.  The 
addition  of  the  notch  (one  more  flaw)  does  not  have  as 
much  impact  on  strength.  But  with  low  a,  the  spread  of 
data  becomes  large.  Designing  to  a  given  degree  of 
reliability,  one  is  forced  to  accept  very  low  stresses. 
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This  is  due  to  the  equation  relating  reliability 
"R" (a)  to  a  and  B  for  an  unnotched  specimen. 

In  a  =  I  ln[-ln  R(a)]  +  In  6  1.24 

where  a  is  the  design  allowable  for  the  desired 
reliability  R.  (see  Appendix  A  Part  II)  For  example,  for 

6=  5,000 

R  =  .9999  (Probability  of  failure  =  1/10,000) 


a 

a  allowable 

a  (average) 

10 

1990 

4760 

30 

3680 

4920 

The  low  a  material  can  use  only  a  small  percentage  of  its 
average  unnotched  failure  stress. 

IV  The  assumed  form  1.13  for  the  variation  of  ineffective 
length  with  notch  size  used  in  the  point-stress  failure  model 
is  confirmed  by  the  Weibull  integrals.  The  form  was  a 
matter  of  convenient  assumption,  and  is  now  satisfactorily 
proven  by  the  integral  results.  See  Tables  1  and  2. 

V  Future  Work  Areas 

The  data  base  for  notch  sensitivity  constants  k  and  m 
and  tensile  data  scatter  a  needs  to  be  expanded.  This  data 
allows  a  fuller  exploration  of  the  brittle  failure  theory, 
and  also  serves  the  designer. 
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The  material  constants  k,  m  need  to  be  related  for 
various  notch  geometries.  This  could  be  obtained  through 
a  study  of  the  stress  solutions  for  general  notch  shapes, 
and  these  solutions  decayed  to  the  specific  geometries 
of  interest-hole,  crack,  etc. 
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Part  2  Biaxial  Tension 


2.1  Weibull  Theory 

A  numerical  integration  can  be  used  to  predict 
basically  conservative  notched  strengths  for  a  material 
subjected  to  biaxial  tension.  The  geometry  studied  was 
that  of  a  circular  hole  in  an  infinite  plate,  and  the 
material  was  assumed  isotropic. 


The  direction  of  failure  cannot  be  presumed  as 
it  was  in  the  case  of  uniaxial  tension.  We  assume  that 
the  most  important  stress  is  the  opening  (mode  1)  stress 
Oq  .  See  Figure  6.  Under  application  of  stresses  0^  , 
far  away  from  the  hole,*  the  tangential,  normal,  and  shear 
stress  shown  in  Fig.  6  are  given  by  [Savin] 

0  =  [(1-p^)  +  (l-4p^+3p^)  cos  20]  2.1 

P  ^ 

+  ^  [(1-p^)  -  (l-4p^+3p‘^)  cos  20] 

*^x  2  4 

00  =  -^  [(1+p  )  -  (l+3p^)  cos  20  2.2 

+  ^  [(1+p^)  +  (l+3p^)  cos  20] 


*  Should  the  original  problem  have  a  non-zero 
this  may  be  reduced  to  a  biaxial  stress  state 


through  an  angle  S,  where 


3  =  ^  tan~^ 


'  2t  ' 

xy 

0 '  -0 ' 

.  X  y 


t'  component, 
by^rotating 


The  "prime"  quantities  are  those  of  the  original  problem. 


FIGURE  6 


STRESS  AND  COORDINATE  DEFINITIONS 


2.3 


Tpg  =  ~  [l+2p^-3p^]  sin  20  +  [l+2p^-3p'^]  sin  26 

R 

where  p  =  /-x - s- 

/x  +  y 

Note  that  as  the  radius  R  0 ,  all  p  terms  vanish,  and 
we  are  left  with  the  standard  stress  transformation 
(rotation  of  axis)  formulae.  [Crandall  &  Dahl] 


The  Weibull  integrals  used  will  be  for  the  opening 
mode  stress  a„  .  For  a  reliability  function  as  in 

— —  O 

uniaxial  tension-  (See  equation  1.5,  1.6) 


where 


B(a) 


/ 

V 


dV 


2.5 


We  again  choose  a  high  risk  volume;  for  biaxial  tension, 
however,  we  must  investigate  a  number  of  rotation  angles  6. 
The  volume  has  length  L,  width  6,  and  plate  thickness  h. 

See  Figure  7.  The  length  L  is  chosen  to  be  long  enough 
to  avoid  finite  width  correction  effects.  The  width  6 
is  chosen  "small"  so  that  the  stress  is  constant  throughout 
the  integration  volume . 


We  assume  that  X  is  the  principal  loading  axis, 
and  divide  equation  2 . 2  through  by  a  : 

A 
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2,6 


^  (1+P^)  +  (l+3p^)  cos  20 


a 

where  r  =  -^  =  LOADING  RATIO 


Evaluating  equation  2.5  for  B(a),  we  use  equation 
2.6,  and  the  fact  that  the  volume  element  dV  =  6h  di^  where 

C 


/~2~2 
=  /x  +y. 


B  (  0  ) 

X 


L 

h6  / 
R 


0 


2.7a 


^  T+l  R  ^  r"  — 1  R  ^  n 

=  h6  ^  (l+(^)  )  +  ^(1+3  (^)  )cos  20]°‘dc 


nt  '  2 

S  R 


2.7b 


We  combine  the  integral  results  and  the  material  shape 
constant  S  into  a  single  shape  constant  for  notched  biaxial 
failure 


This  is  useful,  as  we  can  relate  the  scale  parameter  to 
the  average  notched  failure  strength  by  using  (See  Appendix 
A. I) 

^BF  ""  ^BF  2.9 

Note  that  equation  2 . 6  has  inherently  unknown  constants 
such  as  6  within  it;  by  normalizing  the  notched  strength 
to  an  equivalent  volume  unnotched  strength,  we  can  get 
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rid  of  the  unknown  constants.  Since  we  do  have  an 
isotropic  material,  we  can  compute  an  unnotched  biaxial 
strength  based  also  on  stresses.  Simply  let  p  ->■  0 
in  equation  2.8,  (no  hole) 


Bg  =  6(h6L[^^  +  cos  29]“ 


-1/a 


2.10 


Again,  the  scale  parameter  is  related  to  the  average 
unnotched  strength  by 


Fg  =  8g  r (1+1/a) 


2.11 


where  a_  =  the  average  value  of  a  on  failure . 

15  X 


We  now  ratio  unnotched  strength  to  notched 
strength  to  eliminate  some  of  the  unknown  constants: 


^BN 

“b 


^  y+1  R  ^  T*— 1  R  ^ 

/  [^(l+(^)  )  +  i^(l+3(|)  )cos  20]  dr. 

T-l/arr+1  ,  r-1  ooi"! 

L  '  [— o—  +  — ^  cos  26] 


ri/ct 


2.12 


NOTE  that  this  assumes  that  a  stays  constant  for  the 
unnotched  and  notched  strength  tests.  All  items  in 
equation  2.12  may  now  be  calculated  for  various  angles  f 
at  a  given  loading  ratio  r.  The  angle  0  which  provides 


the  lowest 


and  value. 


BN 

“b 


indicates  the  failure  location  (direction) 
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PLATE  THICKNESS  *  h 


FIGURE  7 

VOLUME  OF  INTEGRATION  FOR  NUMERICAL  INTEGRAL 
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2.2  Numerical  Integral  Results  (Circular  hole) 

The  numerical  integration  necessary  to  evaluate 
the  average  notched  strength  was  performed  by  the  IMSL 
routine  DCADRE.  The  numerical  values  used  in  this  study 
were : 

hole  radius  R  =  0.5" 

integration  length  L  =  5 . " 

Weibull  shape  parameter  a  =  10. 
a 

The  loading  ratio  r  =  —  was  varied  from  0  to  10.  The 
result  for  r  =  0  =  correctly  predicts  the  uniaxial 
strength. 

Table  3 


r 

®fail 

0 

VO 

o 

• 

0.568 

0.1 

90. 

0.584 

in 

o 

90. 

0.659 

1. 

ANY 

0.764 

2. 

0. 

0.659 

3. 

0. 

0.626 

5. 

0. 

0.602 

10. 

0. 

0.584 

The  failure  direction  is  thus  seen  to  "snap 
through"  so  that  it  is  always  perpendicular  to  the 
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largest  applied  load.  This  helps  avoid  awkwardness  in 
the  nature  of  cr^  ,  the  unnotched  strength;  if  6  were 
not  0°  or  90°,  there  is  no  assurance  that  the  unnotched 
sample  would  ever  fail  in  the  direction  9.  Thus,  biaxial 
unnotched  strength  is  equivalent  to  the  uniaxial 
unnotched  strength  "a^"  .  At  r  =  1  we  reach  a  state  of 
equal  stresses  for  any  9 ,  so  the  failure  may  propagate 
in  any  direction. 

The  addition  of  the  second  load  axis  is  also 
seen  to  contribute  to  strength.  This  is  explained  as 
follows.  With  a  single  axis  load,  say  a  only,  there 
are  high  tensile  stresses  at  9  =  90°.  The  addition 
of  Cy  provides  compressive  stresses  which  partially 
offset  the  high  tensile  stresses  in  this  critical  region 
at  9  =  90°. 

The  numerical  integrals  were  performed  for  a 
series  of  hole  sizes  and  compared  with  data  from  a 
quasi-isotropic  graphite/epoxy  laminate  [0/±45/90]g 
[Daniel] .  The  results  are  shown  first  for  uniaxial, 
then  biaxial  strengths.  The  a  value  is  not  known,  but 
is  probably  in  the  20-30  range  from  experience. 
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Table  4 


Hole  Diameter 
(2R)  Inch 

Exper.  Notched  Str. 
CTjj  ,  MPa 

Exper . 
<^n/'^o 

Calculated 

a  =  10 

cr„/a 
N'  o 

a  =  25 

1.0 

219 

0.436 

0.568 

0.429 

0.75 

231 

0.460 

0.584 

0.434 

0.50 

255 

0.508 

0.608 

0.441 

0.25 

276 

0.550 

0.651 

0.454 

The  integrals  for  the  calculated  used  integration 

L  =  5",  loading  ratio  r  =  0.  No  finite  width  correction 
factors  were  used  in  experimental  or  calculated  results. 
The  experimental  results  are  bracketed  by  calculated 
values  for  a  =  lO  and  a  =  25.  The  caculations  for  a  =  25 
are  conservative. 


Table  5 


Hole  Diameter 

2R  Inch 

Exper.  Biax.  Str. 
<7bf  ^  MPa 

Exper. 

’bf/Ob 

Calculated 

a  =  10 

a  =  25 

1.0 

276 

0.550 

0.764 

0.621 

0.75 

280 

0.558 

0.785 

0.629 

0.50 

320 

0.637 

0.814 

0.639 

0.25 

366 

0.729 

0.862 

0.657 

We  see  for  the  first  time  a  non-conservative  failure 
prediction  Cgp/Cg  with  the  calculation  up  to  13%  higher 
than  experiment.  However,  the  data  scatter  in  the 
experimental  results  is  large,  and  the  results  somewhat 
inconclusive.  (Experimental  result  =  0.558  is  average 

of  two  points  0.614,  0.543.) 
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2.3  Conclusions  and  Recommendations 


I  The  calculated  biaxial  failure  loads  are  generally 
conservative  versus  experimental  data.  That  is,  the 
numerical  Weibull  integrals  predict  a  more  severe 
strength  degradation  in  a  notched  sample  than  we  actually 
experience.  The  few  cases  of  non-conservatism  (see  Table 
5)  are  inconclusive  due  to  a  scarcity  of  data.  The 
conservatism  of  calculation  occurs  as  the  calculations 
use  the  highest  risk  volume  with  a  "weak  link"  theory. 

In  actuality,  there  would  be  some  "give"  to  the  material. 

II  The  Weibull  biaxial  integration  results  accurately 
predict  the  strengthening  effect  of  the  stress  along 

the  second  axis.  This  strengthening  effect  is  due  to 
compressive  stresses  from  one  axis  loading  partially 
cancelling  the  tensile  stresses  from  the  other  axis 
loading . 

Further,  the  Weibull  integration  predicts  within 


15%  the  magnitude  of  the  strengthing 

effect 

for  equal 

axis  loading: 

Hole  Diameter 
2R  Inch 

(from  Tables  4,  5) 

Experimental 

(a  =a  ) 

X  y 

Calculated  a„_/G 

or 

(a^=Oy) 

a-10  0=25 

1.00 

1.26 

1.35 

1.45 

0.75 

1.21 

1.34 

1.45 

0.50 

1.25 

1.34 

1.45 

0.25 

1.33 

1.32 

1.45 

AVG 

1.26 

1.34 

1.45 
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where  a__/a  is  the  biaxial  strength  with  a  =  a 
BF'  o  ^  X  y 

divided  by  the  uniaxial  strength. 

III  The  Weibull  integration  results  for  the  failure 
direction  agree  with  the  limited  amount  of  experimental 
data.  Experimentally,  the  biaxial  samples  filed  perpendi¬ 
cular  to  the  stronger  load,  with  inconclusive  results 

for  the  equal  biaxial  loading  case.  [Daniel]  Some  of 

the  materials  tested  in  reference  [Daniel]  were  anisotropic, 

and  no  predictive  calculations  have  been  performed. 

IV  Future  Work  Areas ; 

Fracture  of  anisotropic  materials  in  biaxial 
tension  can  be  similarly  addressed  by  numerical  evaluation 
of  Weibull  integrals.  The  failure  direction  of  such 
materials  can  be  predicted,  along  with  some  measure  of 
strength  reduction  in  the  presence  of  a  notch.  These 
calculations  should  be  supported  by  more  biaxial  failure 
data.  Having  done  this  for  the  biaxial  load  case,  it 
should  be  simple  to  do  the  uniaxial  load  case  as  a 
reduction. 

Experience  needs  to  be  gained  in  how  to  properly 
test  samples  in  biaxial  loading.  The  simple  tube 
pressurization  test  is,  of  course,  not  suitable  for 
notched  geometries.  The  existing  literature  dealing 
with  this  problem  is  very  small . 
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APPENDIX  A:  Mathematical  Development  of  the  Weibull 
Distribution 

I  Mean  and  Standard  Deviation 

The  generally  used  form  for  a  cumulative  distribution 

function  (or  simply  "distribution  function")  F  is: 

F(a)  =  PCCg  _<  a)  I.l 

where  a  is  a  random  variable*,  P  denotes  probability,  and 
s 

a  is  the  particular  value  of  interest.  For  a  set  of  data 
which  obeys  a  Weibull  distribution,  the  distribution 
function  is: 

F(a)  =  1  -  exp [- (a/3) °‘']  a  ^  0  1.2 

where  exp  is  exponential  function 

3  is  the  scale  or  location  parameter 
a  is  the  shape  parameter 

Mathematicians  often  use  the  equivalent  notation 

a  ~  W(3,a)  1.3 

to  show  that  a  obeys  a  Weibull  distribution. 

We  are  often  interested  in  the  reliability,  i.e. 
the  probability  that  the  random  variable  (outcome  of  an 
experiment)  exceeds  the  particular  value  of  interest. 

This  reliability  is  expressed  as 


*A  "random  variable"  is  simply  the  outcome  of  an  experiment. 
For  example,  the  value  shown  on  a  fair  die  after  a  die  toss 
is  a  random  variable. 
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1.4 


R(a)  =  P(a  >  g) 

S 

or  R(a)=l-F(a) 

and,  for  a  Weibull  distribution,  is  given  by; 

R(a)  =  exp  [- (0/3)  0^0  1.5 

We  now  proceed  to  find  the  mean  or  average  in  terms 
of  the  parameters  a, 3.  To  get  the  probability  density 
function  (p.d.f.)  p(0)  of  a  distribution,  take  the 
derivative  as  follows: 

Since  p.d.f.  p{0)d0  =  P[0-do  <0  <  cr  +  d0]  1.6 

P(cr)  =  ^  tF(0)]  1.7 

Using  the  Weibull  form  for  the  distribution  function  F(I.l), 

p(0)  =  3‘“  exp[-(0/3)°‘]  1.8 

The  formula  for  obtaining  the  mean  is : 

00 

E[0_]  =  0  =  /  0  p(0)  d0  1.9 

o 

where  F  denotes  the  average,  and  E  is  the  expected  value 
operator.  The  form  for  the  Weibull  p.d.f.  1.8  is 
substituted  into  1.9.  The  ensuing  integration  uses  the 
change  of  variable 

o“06  a 

X  =  3  cr 

which  does  not  change  the  limits  of  integration.  The 
final  expression  is 


34 


a  =  ^  f  e“^  dx  1. 10 

o 

This  form  is  identical  to  the  Gamma  Function  T,  a  table 
look-up  in  mathematical  handbooks .  Thus  we  have 

a  =  er (1+1/a)  I. 11 

For  a^l,  O.SSe^r^l.  ,  so  that  o’  £  B;  or  the  average 
is  less  than  the  scale  parameter  due  to  the  skewness  of 
the  distribution. 

The  standard  deviation  for  a  Weibull  distribution 
may  be  found  similarly  from  the  basic  definition: 

St. Dev.  S  =  {E(a^]  -  (E[a^])^}^  1.12 

s  s 

When  the  definition  1.12  is  used  with  the  p.d.f.  definition 
1.8,  it  yields: 

S  =  {3^  [r (1+2/a)  -  (1+1/a) 1.13 
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II  Statistical  Inference  of  Weibull  Parameters 

Given  a  set  of  data  from  an  experiment,  we  would 
like  to  predict  point  extimates  a  and  d  for  the  entire 
(infinite)  population.  One  way  to  do  this  is  by  use  of 
logarithms;  this  method  should  be  used  only  if  a  digital 
computer  or  programmable  calculator  is  not  available. 

The  second  method  is  that  of  solution  of  the  maximum 
likelihood  equations;  this  demands  use  of  a  computer 
(or  programmable  calculator) . 
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We  proceed  to  manipulate  the  formula  1.2  into  a 
convenient  linear  form.  Since  P(a  <  a)  =  F(a), 

S 

1-F  =  exp[-(a/6)°‘]  II. 2 

Taking  natural  logarithms, 

In(l-F)  =  -(a/3)°‘  II. 3 

Clearing  the  minus  sign  and  again  taking  natural 
logarithms , 

ln{-ln(l-F)  =  aln(a/3)  II. 4a 

or 

Ina  =  i  ln[-ln(l-F)]  +  ln3  II. 4b 

The  equation  II. 4b  is  in  linear  form  for  a  linear  squares 
analysis  of  data  points  Ina  versus  In [-In (1-Rank) ] .  The 
slope  is  1/a,  the  intercept  is  ln3. 

Note  that  the  form  II. 4b  is  also  usable  to  calculate 
an  allowable  a  for  a  given  structural  reliability  R  =  1-f, 
once  the  a  and  3  are  known. 

Problems  with  logarithm  method ;  The  a , 3  found  by 
linear  regression  on  II. 4b  are  point  estimates.  To  obtain 
meaningful  interval  estimates  of  say  a  requires  us  to 
assume  that  the  least  squares  statistics  1/a  and  Ind  are 
normally  distributed.  This  assumption  seems  doubtful  at 
best.  The  logarithm  method  also  has  difficulty  handling 
censored  data  points,  which  can  occur  in  fatigue  experiments 
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(as  run-out,  etc.)  . 


Maximum  Likelihood  Method 

This  method  has  a  sound  theoretical  basis  and  allows 
interval  estimates  to  be  made  for  a,d.  However,  it 
demands  solution  of  a  very  non-linear  equation. 

/N 

The  maximum  likelihood  estimator  6  for  a  set  of 

n 

n  random  variables  is  given  by  the  solution  of : 

L(0)  =  f(x  ;9)f(x^;e)  ...  f(x„;6)  II. 5 

•  ^  Xi 


— 

n 

— 

fX  .  1 

a-1 

rX  . 

a 

MAX  L(0)==  _ 

J 

ot 

-  3 

1 

J 

exp 

— 

.± 

II. 6 

where  f  is  the  p.d.f.  of  X,  and  9  is  the  true  population 

/V 

arameter,  and  L(6)  reaches  a  maximum  at  9  .[A.l]  maximi¬ 
zation  equation  II. 6  acts  to  maximize  the  likelihood  of 
"legitimizing"  the  outcome  of  a  given  experiment . [A. 2] 


The  maximum  likelihood  equations  for  d,3  are  given  by: 


n 


k(d)  =  0  = 


I  X“  InX. 
i=l  ^  ^ 


n 


l  InX. 
1  i=l  ^ 


.1  ^i 


S 


n 


II. 7 


i=l 


3  = 


,  n  ^ 
i  I  X? 
"  ill  ^ 


1/a 


II. 8 


/V 

The  d,6  are  point  estimates  to  the  true  population  parameters 
a, 3.  [A. 3]  In  practice,  equation  II. 7  is  solved  for  d  and 
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the  result  used  in  II. 8  to  get  3.  For  small  sample  sizes, 

/s. 

the  6t,3  estimates  from  equations  II. 7,  II.  8  are  biased. 
Fortunately,  there  are  correction  factors  available  which 
depend  only  on  n,  the  number  of  data  points  [A. 4].  These 
correction  factors  B(n),  shown  below,  thus  provide  unbiasing 
so  that 

E  [  B(n)  a  ]  =  a  II. 9 


TABLE  A.l 

Unbiasing  Factors  B (n)  for  MLE  of  a 


n 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

B(n) 

.669 

.752 

.792 

.820 

.842 

.859 

,872 

.883 

.893 

.901 

.908 

.914 

n 

18 

20 

22 

24 

26 

28 

30 

32 

34 

36 

38 

40 

B(n) 

-923 

.931 

.938 

.943 

,947 

.951 

.955 

.958 

.960' 

-962 

.964 

.966 

n 

42 

44 

46 

48 

50 

52 

54 

56 

58 

60 

62 

64 

B{n) 

.968 

.970 

.971 

.972 

.973 

.974 

.975 

.976 

.977 

.978 

.979 

.980 

n 

66 

68 

70 

72 

74 

76 

78 

80 

85 

90 

100 

120 

B(n) 

.980 

.931 

.981 

.982 

.982 

.983 

.983 

.984 

.985 

.986 

.987 

.990 
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Interval  estimates  for  a  and  B  may  be  obtained 
for  any  desired  level  y  of  significance  desired.  [A. 4] 
The  intervals  can  be  one-sided  or  two-sided. 
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APPENDIX  B:  Physical  Development  of  the  Weibull  Distribution 


The  Weibull  distribution  was  developed  as  a  distri¬ 
bution  of  wide  utility  for  fitting  data.  It  has  a  parti¬ 
cular  physical  significance  in  treating  fracture  data  of 
brittle  material.  The  central  thesis  of  this  concept  is 
that  the  strength  of  a  material  cannot  be  expressed  by 
a  single  number.  The  risk  of  breakage  may  vary  throughout 
the  material;  the  break  site  need  not  be  the  site  at  which 
the  ultimate  stress  is  first  reached. 

If  the  material  is  viewed  as  a  series  of  n  volume 

elements  and  the  individual  element  failure  probabilities 

are  S, ,  S„,...S  ,  then  the  survival  of  the  whole  (no  failure) 
1  2'  n 

depends  on  the  survival  of  all  n  elements,  i.e. 

1-S  =  (1-S3^)(1-S2)  ...  (1-S^)  1.1 

NOTE  the  "weak  link"  effect  implied  here;  failure  of  any 

subpart  brings  failure  of  the  whole.  For  infinitely  small 

elements,  the  probability  of  failure  in  any  one  given 

element  becomes  vanishingly  small,  thus,  in  1.1  we  multiply 

the  right-hand- side  and  omit  the  higher  order  terms: 

S  =  S,S-S,  ...  1.2 

1  2  o  n 

The  probability  of  an  individual  volume  failure  is 
proportional  to  the  volume;  letting  F(x)  be  the  distribution 
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function  for  the  infinitely  small  volume. 

Si  =  F^(a)  dV  i  =  1,2, ...n  1.3 

where  is  finite  and  dV  is  infinitely  small. 

The  value  of  a,  which  could  represent  stress  or  strain,  is 
in  general  a  functicn  of  position,  i.e.  a  =  a{X).  If 
*^MAX  the  maximum  stress  at  any  point  in  the 

material,  then  the  local  stress  may  be  considered  a 


product  of  and  a  geometric  locus  function  t(x)  so  that 


"  *  “max  - 


From  1.1,  we  obtain  generally  that 


n 

log(l-S)  =  I  loq(l-S)  1.5 

M=1 

As  the  number  of  volume  elements,  n,  increases  indefinitely, 
converges  to  zero,  so  that 

log(l-Sjj)  =  1.6 

Using  equation  1.5  as  n  ^  , 


log (1-S) 


-lim 

n^oo 


n 


1.7 


In  the  limit,  using  1.3  with  1.7, 


log(l-S)  =  -/  F(a)dV 
V 

If  we  define  for  simplicity 

B  =  -/  F(a)dV 
V 


1.8 


1.9 
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Then  the  distribution  function  for  the  entire  body 

F(a)  =  1  -  1.10 

This  F  gives  the  probability  of  failure  for  the  outcome  of 
any  particular  experiment  as  a  function  of  the  applied 
stress  a. 


Weibull  [B.l]  chooses  an  elementary  distribution 
function  which  satisfies  equation  1.1 


F(a)  = 


a-a 


u 


1.11 


where  a  is  the  value  below  which  no  volume  will  ever 
u 

fail  (usually  taken  to  be  zero) . 
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